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ABSTRACT 

We apply the image subtraction method to re-analyze the ESO Very Large 
Telescope data on M83 (NGC 5236), obtained and analyzed by Thim et al. 
Whereas Thim et al. found 12 Cepheids with periods between 12-55 days, we find 
112 Cepheids with periods ranging from 7-91 days, as well as ~ 60 other vari- 
ables. These include 2 candidate eclipsing binaries, which, if confirmed, would be 
the first optically discovered outside the Local Group. We thus demonstrate that 
the image subtraction method is much more powerful for detecting variability, 
especially in crowded fields. However, HST observations are necessary to obtain 
a Cepheid period-luminosity distance not dominated by blending and crowding. 
We propose a "hybrid" approach, where numerous Cepheids are discovered and 
characterized using large ground-based telescopes and then followed-up with the 
HST to obtain precise distances. 

Subject headings: Cepheids — distance scale — galaxies: individual (M83, NGC 
5236) 



1. Introduction 



Cepheid variables are one of the most important primary distance indicators due to 
the period-luminosity (PL) relation they obey. Precise Cepheid distances require a large 
number of Cepheids with accurate photometry and periods. However, the more distant 
the galaxy, the more important become the effects of crowding and blending with nearby 
luminous stars (Mochejska et al. 2000; Stanek & Udalski 1999; Mochejska et al. 2001), which 
introduce systematic, one-sided errors in the distance measurements and skew the distance 
towards artificially lower values. These effects can only be disentangled satisfactorily with 
the Hubble Space Telescope, which has the required spatial resolution, especially now with 
the new instrumentation (ACS). 
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The spiral galaxy M83, the principal member of the nearby M83 group, is a good 
target to obtain a Cepheid distance to. Observations of M83 with the Very Large Telescope 
(VLT) were analyzed by Thim et al. (2003) with point-spread function (PSF) photometry, 
discovering 12 Cepheids, which were used to derive the Cepheid distance to the galaxy. 
However, the "traditional" method of doing PSF fitting photometry on all the stars on an 
image and looking for variations in the light from night to night is not very efficient or 
effective in finding Cepheids in crowded fields. In this paper we re-analyze the VLT data for 
M83 using image subtraction to demonstrate this fact. 

Observations of M83, at a distance of ~ 4.5 Mpc (Thim et al. 2003), with the 8.2 meter 
VLT and 0.76" median seeing are roughly equivalent to the DIRECT project observations 
of M31/M33 with the FLWO 1.2 meter telescope, at a distance of ~ 780 kpc (Stanek & 
Garnavich 1998) with ~ 4" seeing. This fact partly motivated this paper - we would con- 
sider such poor seeing data unacceptable to run PSF photometry on. The three most recent 
DIRECT papers (Mochejska et al. 2001a,b; Bonanos et al. 2003) have used Alard's image 
subtraction package ISIS (Alard & Lupton 1998; Alard 2000) to discover variables, in partic- 
ular detached eclipsing binaries and Cepheids, in M31 and M33. This method has become a 
method of choice for variability searches in crowded fields. 

At the distance of M83, the issue of blending must be taken into account in deriving the 
Cepheid distance. The median seeing of the VLT data is 0.76", which corresponds to 17 pc 
in M83. As first discussed by Mochejska et al. (2000), blending is the close association of 
a Cepheid with one or more intrinsically luminous stars, which is the result of the higher 
value of the star-star correlation function for massive stars, such as Cepheids, compared to 
random field stars. This effect cannot be detected within the observed PSF by usual analysis. 
Mochejska et al. (2000, 2001c) address the effect of blending on the Cepheid distances to M31 
and M33. They compare high resolution HST images to the ground-based DIRECT data 
and find that blending can affect the flux of a Cepheid typically by ~ 20 — 30%, which leads 
to an underestimation of the true distance. This phenomenon is different from crowding, 
which is the random background luminosity fluctuation in each resolution element. In M83, 
a large fraction of the flux of a blended Cepheid could come from its companions, assuming 
the effect scales to that in M31 and M33. Based on a simulation by Stanek & Udalski (1999), 
this would result in a significant distance bias. Whereas the discovery of Cepheids in nearby 
galaxies can be done adequately from the ground given good signal-to-noise photometry, 
deriving the Cepheid PL distance requires high spatial resolution HST imaging. 

In this paper, we apply the image subtraction method to the M83 dataset. In Section 
2 we describe the observations, in Section 3 the image subtraction method and finally the 
results, which are discussed in Section 4. 
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2. Observations 

We have retrieved the M83 data obtained by FORS1 on the ESO Very Large Telescope, 
from the ESO/ST-ECF Science Archive Facility. The field of view is 6.8 x 6.8 arcminutes, 
with a pixel scale of 0.2" per pixel. The 2048 x 2048 Tektronix CCD has 24 jj,m pixels. There 
are 34 epochs in F-band spanning a period of 1.5 years, from January 2000 to July 2001. 
Each epoch consists of 3-4 subexposures of 400 or 500 sec. More details are given by Thim 
et al. (2003). 

3. The Method of Image Subtraction and Results 

The images were overscan corrected and flat fielded with standard IRAF 1 routines. The 
transformation from rectangular to equatorial coordinates was derived using 84 transforma- 
tion stars from the USNO-B1.0 (Monet et al. 2003) catalog. The average difference between 
the catalog and the computed coordinates for the transformation stars was less than 0."3 
in RA and 0."3 in Dec. Next, we ran the image subtraction package ISIS (Alard & Lupton 
1998; Alard 2000) on the 104 V-band images of M83. 

The ISIS reduction procedure consists of several steps. Initially, all the frames are trans- 
formed to a common coordinate grid. The best seeing frame, FORS.2000-04-05T06:16:02.924 
in our case, is chosen to be the reference image. Next, a composite reference image is cre- 
ated by stacking several best seeing frames. For each frame, the composite reference image 
is convolved with a kernel to match its PSF and then subtracted. On the subtracted images, 
the constant stars cancel out, and only the signal from variable stars remains. A median im- 
age is constructed of all the subtracted images, and the variable stars are identified visually 
as significant peaks. Finally, profile photometry is extracted from the subtracted images. 
Mochejska et al. (2001a) describes this procedure in more detail. 

After locating the positions of the bright peaks, we obtained light curves for ~ 650 
candidate variables. We ran them through the DIRECT pipeline (Kaluzny et al. 1998; 
Stanek et al. 1998) which fits model Cepheid light curves and eclipsing binary light curves 
and classifies the stars as Cepheids, eclipsing binaries (EBs) or other variables. After checking 
these light curves, we found 112 Cepheids, with periods ranging from 7 to 91 days, 2 candidate 
EBs and ~ 60 other variable stars. Figure 1 presents 10 sample Cepheid flux light curves. 
Table 1 lists the coordinates (RA, Dec) and (X, Y) on the reference image for each Cepheid 
and the periods. We recover all but one of the Cepheids found by Thim et al. (2003), 



1 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Asso- 
ciation of Universities for Research in Astronomy, Inc., under cooperative agreement with the NSF. 
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the Cepheid C4. The periods we derive agree well with those of Thim et al. (2003), the 
median difference being 0.04 days and the largest difference being 0.22 days for Cepheid C6. 
Figure 2 plots the location of the Cepheids as circles on the reference image, their size being 
proportional to the period. A 29 day Cepheid is labeled for scale. On the same plot, the 
stars show the position of the 12 Cepheids from Thim et al. (2003). The image subtraction 
method detects many Cepheids in the more crowded part of the field, which is much harder 
to do with the "traditional" PSF photometry method. 

We do not claim completeness in finding Cepheids, however we do show that Cepheids 
in the highly crowded and blended central regions of the galaxy can be detected from the 
ground with high quality data, such as these VLT observations. High resolution observations 
are necessary to disentangle blending for Cepheids both in the outer and inner parts of the 
galaxy. 

4. Discussion 

This paper re-analyzes the excellent VLT data of M83 obtained by Thim et al. (2003) 
using the image subtraction method. The resulting nine-fold increase in the number of 
Cepheids detected indicates that image subtraction should be used in crowded fields. We 
also present parameters and sample light curves for the 112 Cepheids we have found. These 
additional Cepheids are valuable for determining the PL distance to M83 accurately. How- 
ever, HST observations are necessary to resolve blending effects. After we started working 
on this project, we became aware of the Cycle 12 program with the title "M83: Calibrating 
the Cepheid PL Relation" (PI: B. Madore). We note that some of these 112 Cepheids might 
be Population II variables, however, a PL diagram would be necessary to distinguish them. 

In Figure 3, we plot the flux amplitude versus period relation for the 112 Cepheids. We 
have taken the full flux amplitude from the light curve and plotted the logarithm of this 
quantity versus logarithm of the period P. There is a definite correlation, however, it is not 
very tight. Paczyhski & Pindor (2000) demonstrate that the period-flux amplitude relation 
is not universal and needs to be calibrated before being used to measure distances accurately. 

In addition to these Cepheids, we also find ~ 60 other variable stars, including 2 can- 
didate eclipsing binaries (EBs). If these are confirmed from their position on a CMD to be 
located in M83, they would be the first optically discovered EBs outside the Local Group. 
Figure 4 presents some of the more interesting light curves of other periodic or non periodic 
variables that we find, including the candidate EBs. There are two long period variables, 
which are possibly Cepheids with periods of 132 and 195 days. Their positions on the PL 
diagram would verify this. 
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CEPH (13:36:48.0,-29:46:24.9) P=6.777 



CEPH (13:36:50.6,-29:46:40.9) P= 10.587 



3 




0.5 1 1.5 

CEPH (13:36:49.6,-29:49:51.3) P= 12.464 



0.5 1 1.5 

CEPH (13:36:49.3,-29:49:48.6) P = 14.890 




0.5 1 1.5 

CEPH (13:36:52.1,-29:48:56.8) P=15.743 



0.5 1 1.5 

CEPH (13:36:48.2,-29:49:20.7) P= 17.793 




0.5 1 1.5 

CEPH (13:36:53.6,-29:48:52.7) P=20.477 




0.5 1 1.5 

CEPH (13:36:52.9,-29:50:09.5) P=26.343 




0.5 1 1.5 

CEPH (13:36:45.0,-29:44:47.7) P=36.546 



0.5 1 1.5 

CEPH (13:36:4B.l, -29:49:46.8) P=54.831 




Phase 

Fig. 1. — Flux light curves of 10 sample Cepheids found in M83. For each Cepheid, the 
coordinates (RA, Dec) are given in parentheses, along with the period. The median point 
for each epoch is plotted. 
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500 1000 1500 2000 

Fig. 2. — Location of the 112 Cepheids we found with image subtraction on the reference 
image. Circles are centered at the position of the Cepheids and their size is proportional 
to the period of the Cepheid (see labeled 29 day Cepheid for reference). Stars indicate the 
positions of the 12 Cepheids found by Thim et al. (2003). All but one were recovered by our 
analysis. [See the electronic edition of the Journal for a color version of this figure] 
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Fig. 3. — Flux amplitude versus period for 112 Cepheids in M83. 



- 8- 



x 



Var (13:36:44.8, -39:50:01.6) 



Var (13:36:56.6, -29:48:17.5) 




1600 1800 
Var (13:36:56.9, -29:46:20.3) 




1600 1600 
Var (13:36:41.2, -29:48:27.2) 



—I 1 1— 



* » ■ 



_i_ 



— i 1 1 — 



1600 1800 
Var (13:36:47.2, -29:46:41.0) 



1600 1800 

Var (13:36:50.7, -29:45:32.3) 




1600 1800 
Var (13:36:49.7, -29:45:25.8) 




1600 1800 
Var (13:36:42.8, -29:45:16.4) 



T" 



— i 1 1 — 



T 



—i 1 1 — 



i I i i i I i i i 



1600 1800 2000 

Var (13:36:47.7, -29:49:07.8) P= 132.5 



1600 1800 8000 

Var (13:36:53.4, -29:49:23.0) P= 195.1 




0.5 1 1.5 

EB (13:36:57.2, -29:48:16.9) P= 4.95 



-5000 


5000 



1 — T -1 — ' — ' — 1 — T -1 — ' — ' — 1 — T" 1 — r ~ 
. : A" i . i i" i 

j. .■ . ■ ' I. .• 



0.5 1 1.5 

EB (13:36:56.4, -29:47:52.6) P= 9.35 




HJD - 2450000 



Fig. 4. — Flux light curves of some miscellaneous variables in M83, with coordinates (RA, 
Dec). Two possible long period Cepheids and two candidate eclipsing binaries are shown. 
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Table 1. Cepheids in M83 







P 








RA 


Dec 


(days) 


X 


Y 


Comments a 


13:36:50.6 


-29:48:40.9 


10.587 


589.28 


514.04 




13:36:43.4 


-29:49:31.0 


11.313 


1052.49 


264.76 




13:36:42.6 


-29:49:27.0 


11.745 


1106.87 


285.00 




13:36:48.6 


-29:49:31.1 


11.997 


717.45 


263.21 




13:36:49.6 


-29:49:51.3 


12.464 


652.30 


161.81 


C9, 12.47 days 



a Names and periods of Cepheids found by Thim et al. (2003). 



Note. — Table 1 is available in its entirety in the electronic version of 
the Journal. A portion is shown here for guidance regarding its form and 
content. 



